AD-767  186 


STATISTICAL  ANALYSIS  OF  SHORTWAVE 
FADEOUT  OCCURRENCE  FOR  THE  YEARS 
1955  TO  1969 

Richard  G.  Hendl,  et  al 

Air  Force  Cambridge  Research  Laboratories 
L.  G.  Hanscom  Field,  Massachusetts 

28  June  1973 


r 


DISTRIBUTED  BY: 


National  Technical  Information  Service 
U.  S.  DEPARTMENT  OF  COMMERCE 

5285  Port  Royal  Road,  Springfield  Va.  22151 


98  U  91  Of 


AFCRL-TR-73-0385 
28  JUNE  1973 

ENVIRONMENTAL  RESEARCH  PAPERS,  NO.  452 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 

L.  G.  HANSCOM  FIELD,  BEDFORD,  MASSACHUSETTS 


Statistical  Analysis  of  Shortwave  Fadeout 
Occurrence  f  >r  the  Years  1955  to  1969 


RICHARD  G.  HENDL,  Copt,  USAF 
ROBERT  A.  SKRIVANEK 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  SYSTEMS  COMMAND 
United  States  Air  Force 


ACCESSION  lor 


NIIS 

BUG 

RJKiOWCEJ 


iffillllMIMH. 


oisTitisinios/Aiiiuiiun  coses 


Dlsl.  AVAIL.  STKIAL 


Qualified  requestors  may  obtain  additional  , 

Teoh™Tw^rneUce0'hm  Sh°“,d 


copies  from  the  Defense 
I  apply  to  the  National 


Unclassified 
Security  Cl«in[fVc«(ion 


DOCUMENT  CONTROL  DATA  -  RAD 

f Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 


t.  ORIGINATING  ACTIVITY  (Corporate  author) 

Air  Force  Cambridge  Research  Laboratories  (PHS) 
L.G.  Hanscom  Field 

Bedford,  Massachusetts  017  30 

2&  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

20.  "GROUP 

1  REPORT  TITLE 

STATISTICAL  ANALYSIS  OF  SHORTWAVE  FADEOUT 

OCCURRENCE  FOR  THE  YEARS  1955  TO  1969 

4.  oesc  rip  tive  NOTES  ( Type  of  report  and  inclusive  dates ) 


5.  AU T  HO R(5)  (First  nee ne,  middle  initial,  last  name ) 

Richard  G.  Hendl,  Capt,  USAF 
Robert  A.  Skrivanek 


*.  REPORT  DATE 

28  June  1973 

7a  TOTAL 

lb.  NO.  OP  REFS 

2 

6 0>  CONTRACT  OR  GRANT  NO. 

b.  PROJECT,  TASK.  WORK  UNIT  NOS.  86660101 

C.  OOD  ELEMENT  62101F 

d,  OOD  SUBELEMENT  681000 

•  a  ORIGINATOR'S  REPORT  NUMREtfS; 

AFCRL  -TR-73-  0385 

s b.  ot HER  REPORT  H&S)  (Any  other  numbers  that  may  be 

nponT  ERP  No.  452 

10.  DISTRIBUTION  STATEMENT 


Approved  for  public  release;  distribution  unlimited. 

tl  SPONSORING  OILITARY  ACTIVITY  " 

Air  Force  Cambridge  Research 
Laboratories  (PHS) 

L.G.  Hanscom  Field 
Bedford,  Massachusetts  01730 

t3.  ABSTRACT 

Shortwave  fadeout  data  for  the  period  1955  through  1969  were  analyzed 
and  correlated  with  solar  flare  data  for  the  same  period.  Details  of  the 
analysis  are  presented  along  with  various  statistical  presentations  of  the  data. 
The  validity  of  computerised  objective  procedures  for  matching  geophysical 
data  is  shown.  Some  doubt  is  cast  on  the  quality  of  the  shortwave  fadeout  data 
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Abstract 


Shortwave  fadeout  data  for  the  period  1955  through  1969  were  analyzed  and 
correlated  with  solar  flare  data  for  the  same  period.  Details  of  the  analysis  are 
presented  along  with  various  statistical  presentations  of  the  data.  The  validity  of 
computerized  objective  procedures  for  matching  geophysical  data  is  shown.  Some 
doubt  is  cast  on  the  quality  of  the  shortwave  fadeout  data  during  the  early  part  of 
the  period  considered. 
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Statistical  Analysis  of  Shortwave  Fadeout  Occurrence 

for  the  Years  1955  to  1969 


The  intent  of  this  report  is  twofold:  (1)  to  acquaint  interested  researchers 
witli  the  results  of  an  objective  procedure  to  match  shortwave  fadeout  (SWF)  and 
solar  flare  reports,  and  (2)  also  to  summarize  the  statistical  properties  of  the 
SWF  phenomena.  Assuming  anyone  using  the  matched  SWF/solar  flare  data  will 
insist  on  knowing  the  procedure  used  in  matching  the  data,  this  report  contains  a 
discussion  of  the  matching  rules  and  their  underlying  logic. 

Anyone  using  the  parts  of  this  report  that  deal  with  the  statistical  behavior  of 
•SWF  phenomena  should  carefully  consider  the  limitations  of  the  observational  net¬ 
work  and  the  effect  these  limitations  have  on  the  validity  of  the  data  itself.  For 
example,  during  the  period  considered,  the  number  of  observatories  was  increas¬ 
ing  every  year  but  their  geographical  distribution  was  not  uniform.  The  largest 
concentration  occurred  in  Europe  and  North  America.  The  results  of  these  limita¬ 
tions  can  be  seen  in  the  data  reported  here. 
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2.  DATA  SOURCES 

SWF  data  for  the  period  1955-1969  were  obtained  directly  from  the  CRPL-F, 
Part  B  series  (Solar  and  Geophysical  Data),  published  by  the  Institute  for  Tele¬ 
communication  Sciences  and  Aeronomy,  These  listings  contain  composite  events 
constructed  from  concurrent  or  consecutive  observations  made  at  one  or  more  ob¬ 
servation  sites.  The  reporting  stations  are  listed,  and  various  subjectively  deter¬ 
mined  indices  are  provided  to  more  completely  describe  the  reconstructed  event. 
The  "widespread  index"  is  an  indication  of  the  geographical  extent  of  the  SWF, 
while  the  "importance  index"  reflects  a  judgment  as  to  the  severity  of  the  event. 

Prior  to  September  1955,  individual  station  reports  were  published.  Such  re¬ 
ports  (apparently  describing  the  same  event)  were  combined  into  one  composite 
description  in  order  to  i.se  a  consistent  listing.  No  "widespread  index"  or  "impor¬ 
tance  index"  were  assigned  to  these  events. 

Flare  data  used  in  the  study  discussed  in  this  report  consisted  of  a  modified 
version  of  the  flare  list  compiled  by  Fischer  and  Hendl  (1969)  [see  Appendix  A], 

3.  DEVELOPMENT  OF  ASSOCIATION  CRITERIA 

If  the  results  from  previous  studies  are  accepted,  a  SWF  characteristically 
begins  7  min  after  the  flare  begins  but  always  before  the  flare  reaches  its  maxi¬ 
mum  phase  (Smith  and  Smith,  1963).  This  suggests  that  an  adequate  matching  pro¬ 
cedure  need  only  require  that  the  SWF  begin  in  the  interval  between  the  beginning 
and  maximum  phases  of  the  flare.  Unfortunately,  during  early  data  gathering,  fre¬ 
quent  breaks  in  worldwide  flare  patrol  and  the  frequent  absence  of  the  time  of  the 
maximum  phase  of  the  flare  made  such  a  matching  procedure  useful  for  only  a  por¬ 
tion  of  the  cases.  The  procedure  finally  adopted  was  able  to  handle  a  number  of 
difficulties  which  arose  when  analyzing  incomplete  data 

The  first  problem  considered  was  the  determination  of  a  suitable  interval  be¬ 
tween  the  begin  times  of  both  the  SWF  and  the  flare  within  which  a  match  would  be 
made.  The  interval  selected  had  to  be  large  enough  to  include  the  uncertainties  in 
detecting  the  individual  events  but  not  be  so  large  as  to  match  two  unrelated  events. 
Cases  when  the  SWF  was  observed  prior  to  the  onset  of  the  flare  had  to  be  allowed 
due  to  minor  breaks  in  the  flar-’  patrol  and  the  uncertain  detectability  of  flare  be¬ 
ginnings. 

Presumably,  there  are  limiting  time  differences  attained  by  only  a  small  num¬ 
ber  of  physically  paired  flares  and  SWF.  Even  if  these  differences  were  known,  it 
would  be  undesirable  to  utilize  them  as  matching  criteria  since  it  would  increase 
the  chances  of  matching  events  which  were  not  physically  part  of  the  same 
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phenomenon.  The  first  step,  then,  was  to  construct  a  "matching  interval",  which 
would  pair  the  standard  events  with  a  high  degree  of  confidence,  and  devise  another 
procedure  to  handle  the  pairs  with  unusually  large  begin  time  differences. 

The  routine  matching  interval  was  determined  by  first  constructing  a  frequency 
distribution  (Table  1)  of  differences  in  SWF  begin  time  and  flare  begin  time  per 
5-min  interval  (up  to  ±  90  min).  At  large  begin  time  differences  the  number  of 
cases  per  5-min  interval  became  nearly  constant  so  the  matching  interval  was  de¬ 
fined  to  include  only  those  nonconstant  intervals.  As  cun  be  seen  from  Table  I, 


Table  1.  Frequency  of  SWF  and  Flare  Matches 
Per  5-Min  Begin  Time  Difference  Interval 


At  (Minutes) 

<-30 

a  -30 

a  -25 

a  -20 

a  -15 

a  -io 

>  -5 

-0 

No.  Matches 

2. 5  per 

interval 

39 

48 

69 

110 

2011 

386 

910 

At  (Minutes) 

*5 

<  10 

s  15 

<  20 

<25 

£  30 

>  30 

No.  Matches) 

510 

134 

69 

51 

48 

32 

30  per 

interval 

Note:  At  is  positive  (negative)  when  the  flare  begins  earlier  (later)  than  SWF. 


this  required  that  +10  min  >  Flare  Begin  Time  ••  SWF  Begin  Time  >  -  20  min,  or 
the  flare  to  begin  either  no  earlier  than  20  min  prior  to  or  no  later  than  10  min 
after  the  SWF  had  begun.  This  first  criterion  allowed  for  the  minor  breaks  which 
existed  in  the  flare  patrol. 

A  SWF  not  matched  with  a  flare  by  the  above  procedure  could  be  grouped  into 
four  types: 

(1)  a  SWF  beginning  close  in  time  to  a  flare  with  an  inexact  begin  time; 

(2)  a  SWF  beginning  during  a  period  of  flare  patrol  outage: 

(3)  a  SWF  beginning  more  than  20  min  after  the  beginning  of  the  flare  and 
continuing  concurrently  with  the  flare.  This  type  would  be  characteristic  of  flares 
with  a  gradual  rise  to  maximum  where  the  actual  begin  time  cf  the  flare  is  subject 
to  a  large  uncertainty: 

(4)  a  SWF  which  occurred  far  removed  from  any  flare. 


Those  unmatched  SWF's  occurring  relatively  close  to  the  flare  with  inexact 
begin  times  were  given  further  consideration.  The  begin  time  of  the  flare  was  ad¬ 
justed  by  comparing  the  reported  duration  against  the  average  duration  for  flares 
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of  the  same  importance  (Table  2).  When  the  duration  was  shorter  than  average, 
the  difference  was  subtracted  from  the  begin  time  (unless  the  reported  end  time 
was  also  inexact,  in  which  case  one-half  the  difference  was  subtracted).  For 
flares  >  Importance  V  with  longer  than  average  duration,  10  min  were  arbitrarily 
subtracted  from  the  begin  time.  The  pair  was  considered  matched  if  the  adjusted 
time  difference  fell  within  the  matching  interval. 

At  tl'.is  poirf  any  SWF  beginning  when  no  flare  patrol  was  being  conducted  dur¬ 
ing  the  matching  interval  was  removed  from  further  consideration. 


Table  2.  Mean  Flare  Duration  as  a 
Function  of  Flare  Importance  (After 
Jensen  and  Fischer,  1967) 


Flare  Importance 

Mean  Duration 
(Minutes ) 

0 

18 

1 

27 

2 

54 

3 

96 

4 

140 

Those  SWF's  which  occ  r  red  simultaneously  with  a  flare  for  at  least  75  per¬ 
cent  of  the  SWF  durations  w.  re  matched  with  that  particular  flare.  If  more  than 
one  flare  met  this  criterion,  the  one  with  the  greatest  importance  was  rhosen. 

The  75  percent  value  was  selected  after  subjectively  examining  cases  matched 
manually.  Ideally,  the  flare  should  be  in  progress  for  the  entire  SWF  duration, 
but  flares  prematurely  ended  by  patrol  breaks  required  that  a  less  stringent  figure 
be  used. 

The  number  of  matches  made  as  a  result  of  each  criterion  is  shown  in  Table  3. 
The  majority  of  matches  is  the  result  of  considering  the  begin  times  of  both  the 
flare  and  SWF.  While  the  total  number  of  SWF's  matched  as  a  result  of  adjusting 
flare  begin  times  and  checking  for  a  simultaneous  occurrence  is  small,  the  SWF's 
were  significant  since  their  elimination  would  introduce  a  bias  into  the  data  due 
solely  to  the  incompleteness  of  the  flare  reports.  Since  the  .lata  base  being  pro¬ 
duced  will  ultimately  be  used  to  produce  forecast  techniques,  it  was  important  that 
all  possible  matches  were  made. 

The  potential  matches  affected  by  the  incomplete  flare  patrol  were  reduced 
during  the  1965-1969  period  as  evidenced  by  the  fewer  number  of  matches  made 
by  adjusting  flare  begin  times.  Additionally,  unmatched  SWF's  occurring  during 


Table  3.  Summary  of  SWF  and  Flare  Pairs  by  Type  of  Match 


1955-64 

1965-69 

Total  Number  SWF  Matched  with  Flares 

2540 

187  3 

Matching  Interval 

2300 

1712 

Matching  Interval  after  Adjusting 

Flare  Begin  Time 

48 

6 

7  5%  Common  Occurrence 

192 

155 

Number  SWF  Occurring  During 

No  I10(  Flare  Patrol 

164 

19 

Number  SWF  Unmatched 

517 

198 

3221 

207  1 

periods  of  no  flare  patrol  comprised  only  about  10  percent  of  the  total  SWF's  un¬ 
matched  during  the  latter  period,  compared  to  more  than  30  percent  for  the  earlier 
period.  After  considering  the  final  matching  criteria,  which  were  the  simultaneous 
occurrences  of  a  SWF  and  flare,  about  18  percent  of  the  SWF's  still  remained  un¬ 
matched.  Apparently  those  were  the  result  of  a  flare  producing  a  SWF  during  a 
secondary  maximum  phase  or  other  type  of  noncharacteristic  sequence  of  events. 

4  EVALUATION  OF  THE  SWF  AND  FLARE  MATCHING  PROCEDURE 

The  ,->WF  data  published  in  the  CRPL  bulletin*  contained  the  begin  time  of  the 
flare  which  vas  considered  to  have  caused  the  SWF.  The  flare  was  matched  sub¬ 
jectively  and  represented  the  consideration  of  all  available  flare  parameters.  As 
a  check  on  the  validity  of  the  present  procedure,  the  frequency  distribution  of  the 
flare  and  SWF’  begin  time  differences  for  those  SWF  matched  in  both  lists  were 
compared.  The  resulting  histograms  (see  Figure  1)  are  essentially  identical,  the 
major  difference  is  the  approximately  30  percent  in  matches  occurring  in  the 
AFCRIj  list.  This  difference  is  almost  certainly  due  to  the  different  flare  data 
considered.  It  was  concluded  that  the  computerized  matching  techniques  accurately 
reproduced  the  data  as  published. 

If  the  technique  is  judged  by  considering  the  total  number  of  SWF's  paired  with 
flares,  the  results  are  less  encouraging  at  first  glance  since  only  4336  SWF's  out 
of  5249  (82.  5  percent)  are  matched.  However,  upon  closer  inspection,  it  is  appar¬ 
ent  that  most  of  the  unmatched  SWF's  occur  in  the  earlier  years  when  the  flare 
patrol  coverage  was  less  continuous  and  flare  reporting  less  complete. 


*See  Section  2 
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adequately  and  could  easily  be  applied  to  all  future  SWF  and  solar  flare  data.  The 
extensive  manual  handling  of  data  that  is  presently  used  to  produce  matched  SWF 
and  flare  listing  appears  to  be  unnecessary. 

Using  the  final  matched  list  of  SWF  and  solar  flare  occurrence,  several  ob¬ 
servations  were  made  that  disclose  the  presence  of  fundamental  difficulties  within 
the  SWF  observing  system. 

In  Figure  2,  the  total  number  of  SWF's  observed  is  plotted  as  a  function  of  the 
time  of  occurrence  in  UT.  Looking  first  at  the  total  curve  1955-1969,  it  is  imme¬ 
diately  apparent  that  the  number  of  SWF's  that  were  observed,  peak  at  approxi¬ 
mately  17  00  UT.  There  is  no  geophysical  reason  why  more  SWF's  should  occur 
at  one  time  of  day  than  another.  The  obvious  explanation  for  this  imbalance  is 
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Figure  2.  Total  Number  of  Shortwave  Fadeouts  Observed 
as  a  Function  of  Time  of  Day  (UT) 
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found  in  the  geographical  distribution  of  SWF  observatories;  the  largest  number 
being  in  North  America  and  Europe.  As  more  SWF  stations  were  added  and  more 
emphasis  was  placed  upon  gathering  "gocd"  data,  the  peak  at  17  00  hours  became 
less  pronounced.  This  can  be  seen  from  che  two  other  curves  on  Figure  2,  show¬ 
ing  1955-1964  data  and  1965-1969  data.  Figure  3,  a  plot  of  matched  SWF  versus 
time  of  occurrence,  shows  similar  characteristics. 
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Figure  3.  Number  of  Shortwave  Fadeouts  Matched  with  Solar  Flares 
as  a  Inunction  of  Time  of  Day  (UT) 


Another  example  of  the  improvement  within  the  observation  network  in  recent 
years  is  shown  in  Figure  4.  Here,  the  curves  showing  total  SWF's  and  matched 
SWF's  are  shown  for  the  years  1955  through  1969.  Note  first  that  the  period  of 
the  curve  correlates  well  with  the  Zurich  sunspot  number,  but  that  the  total  num¬ 
ber  of  SWF's  in  the  1969  solar  maximum  was  very  nearly  the  number  observed 
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Figure  4.  Total  Number  of  Shortwave  Fadeouts  and  Number  of  Matched 
Shortwave  Fadeouts  Plotted  for  the  Years  1955-1969.  For  comparison, 
the  Zurich  sunspot  number  if  also  shown 


during  the  1958  solar  maximum  in  spite  of  the  fact  that  the  level  of  solar  activity 
during  the  1969  peak  was  greatly  reduced  from  the  1958  peak.  From  Figure  5, 
in  which  the  total  number  of  Importance  1,  2,  and  3  SWF's  is  plotted  against  time 
of  occurrence  (UT),  another  example  of  bias  introduced  into  the  data  can  be  seen. 
Note  that  the  distribution  of  SWF's  for  the  larger  events  (Importance  3)  is  rather 
flat  — indicating  that,  if  the  SWF  is  sufficiently  large,  most  stations  will  detect  it. 
Smaller  SWF's  (Importance  1)  are  more  easily  missed  and,  hence,  the  detectabil¬ 
ity  is  largely  a  function  of  the  number  of  stations  looking  for  the  phenomena. 

Another  parameter  used  in  describing  SWF  is  the  Wide  Spread  Index  (WSI) 
which  ranges  from  1  to  5.  This  index  is  a  subjective  indication  of  both  the 


geographical  extent  over  which  the  SWF  was  observed  and  the  degree  of  confidence 
placed  in  the  identification  of  the  event.  Most  SWF's  have  a  WSI  of  either  1  or  5 
(Figure  6)  which  leads  one  to  question  the  usefulness  of  this  index. 

It  is  interesting  to  note  that  in  the  1955-1956  data,  74  percent  of  the  SWF  with 
a  SWI  of  5  were  Importance  3  or  greater,  while  in  1966-1967  dat  <  the  number 
dropped  to  only  35  percent.  This  again  indicates  r.n  upgrading  of  the  network, 
more  sensitive  equipment,  more  experience  in  picking  out  the  SWF's  by  the  ob¬ 
servers,  and  consequently  attributing  a  higher  degree  of  confidence  in  their  obser¬ 
vation  and  analysis.  The  distribution  of  SWF's  by  importance  is  further  illustrated 
in  Figure  7. 
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Figure  6.  Distribution  of  Short¬ 
wave  Fadeouts  as  a  Function  of 
Their  Reported  Widespread  Index 


Figure  7.  Distribution  of  Short¬ 
wave  Fadeouts  as  a  Function  of 
Their  Reported  Importance 


Finally,  in  Figure  8,  the  number  of  SWF's  which  occurred  are  plotted  as  a 
function  of  month  of  occurrence.  The  curves  show  still  another  curious  anomaly 
in  which  the  greatest  number  of  SWF's  occurred  near  the  spring  and  fall  equinoxes. 

In  summary  then,  it  should  be  repeated  that  the  SWF  data  from  the  years 
1955-1969  are  inconsistent,  and  consequently  any  conclusions  reached  should  be 
very  carefully  considered. 
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Figure  8.  Total  Nimber  of  Shortwave  Fadeouts  Plotted 
ap  i  Function  jf  lVlonHi  of  Occurrence 


The  1DG5-J969  data  appears  to  he  substantially  more  consistent  than  the 
earlier  data,  bid  it  is  not  possible  to  determine  its  accuracy.  It  is  unfortunate, 
however,  that  the  SWF  observing  network  has  not  thought  it  appropriate  to  intro¬ 
duce  more  modern,  objective,  and  quantitative  means  of  detection  and  measure¬ 
ment  or  these  very  important  phenomena. 
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Appendix  A 


Details  of  Flare  List  Modifications 

The  method  devised  to  construct  flare  events  by  computer  (Fischer  and  llendl, 
'..-JOG)  was  modified  to  produce  a  better  estimate  of  the  event  end  time. 

Originally  the  latest  end  time  appearing  in  the  g.'oimed  flare  reports  was  cho¬ 
sen  as  the  event  end  time.  This  practice  was  responsible  for  producing  a  sizeable 
percentage  of  the  cases  which  disagreed  with  the  Me  Math- 1  lulbert  flare  list  by 
more  than  10  min  (Category  IV  in  the  original  report). 

it  appeared  possihle  to  ohtain  a  better  estimate  of  the  event  end  time  by  a  pro¬ 
cedure  similar  to  that  which  produced  the  event  begin  time.  Unfortunately,  agree¬ 
ment  among  several  observatories  when  dealing  with  flare  end  times  is  poorer  than 
when  considering  begin  times.  This  difference  in  agreement  makes  the  entire  pro¬ 
cedure  impractical.  A  belter  estimate  of  the  event  end  time  was  obtained  bv  con¬ 
sidering  only  the  two  latest  reported  end  times.  If  the  difference  between  the  two 
was  GO  min  or  greater,  the  earlier  times  became  the  event  end  time.  If  the  G0- 
min  interval  was  not  exceeded,  ttie  latest  reported  time  was  used  as  heiore.  The 
reason  for  not  extending  this  practice  to  other  earlier  reported  end  times  is  simi¬ 
lar  to  that  already  described  in  the  discussion  on  the  begin  time  procedure. 

Since  the  logic  which  assigns  the  time  of  the  maximum  phase  of  the  flare  to  the 
event  utilizes  the  end  time,  an  event  witli  a  modified  end  time  might  have  a  differ¬ 
ent  time  of  maximum. 
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